X = linspace(-1, 1);
X = X'
Phi = zeros(size(X), 20);
for i = 1:20
   Phi(:, i) =  X;
end

for j = 0:19
    Phi(:, j+1) = Phi(:, j+1).^j; 
    figure(21)
    plot(X, Phi(:, j+1))
end
hold on;
figure(21)

for k = 1:20
    plot(X, Phi(:, k))
end

hold off;

